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Introduction 



The first documented attempt to construct the geometry theory in an axiomatic way was made, 
as we know, by Euclid (III cent BC) in his Elements. And while the word 'geometry' literally 
means 'earth measuring', Euclidean geometry doesn't describe elliptic space, as more proper 
for measuring of our planet. New axiomatic approach was revolutionary one, however the 
axiomatic has limitations. Euclid study what can be constructed calculated or demonstrated 
starting with compass and straightedge. It was sufficient for that time. However today, despite 
the fact Euclidean geometry is studied in the school, many people, including geometriests, 
can't remember its axioms. Exception makes famous Euclid's V-th postulate, which many of 
us remember in the form: "At most one line can be drawn through any point not on a given 
line parallel to the given line in a plane" . Euclid decided to formulate it so: "If a line segment 
intersects two straight lines forming two interior angles on the same side that sum to less than 
two right angles, then the two lines, if extended indefinitely, meet on that side on which the 
angles sum to less than two right angles" . 

For modern geometry Euclid's axiomatic has several limitations: 

• Euclid's axiomatic theory covered the only geometry system and only two-dimensional 
case. The axiomatic of Euclidian geometry used today was developed by David Hilbert 
(1862 — 1943), has 20 axioms and covers two and three dimensions. 

• The four-dimensional case uses much more axioms. Development of axiomatic for spaces 
of further dimensions is non-trivial. 

• Except hyperbolic geometry, construction of good axiomatic for other geometries is also 
non-trivial. Usually, an axiomatic is constructed after the geometry is well studied with 
aim of some model (for example, [3] describes the space-time axiomatic). 

• Undefined notions in geometry (point, line, between) differ very much from undefined 
notions in other mathematic disciplines (number, function, space). Undefined notions of 
different geometries differ from each other. 

• Mathematicians successful study Euclidean space of any dimension using analytic geom- 
etry and forget Euclid's axiomatic. 

Euclid's axiomatic played one important role. Its V-th postulate is so hard expressed and 
creates so artificial feeling that urged mathematicians to create the hyperbolic geometry. Sad, 
when Nikolai Lobachevsky (1792 — 1856) and Janos Bolyai (1802 — 1860) published their 
results, the new geometry was slow in acceptance. Only after decades it was demonstrated 
that hyperbolic geometry is interior geometry of surfaces with constant negative curvature. 
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Introduction 



After next several years some models of hyperbolic geometry were elaborated. Due to that 
fact the new geometry became accessible. 

Author of a model, Felix Klein (1849 — 1925) proposed "Erlangen Program" [ | - the 
unified view over different geometries as complex of different transformation groups of space. 
The invariants of these groups are figures of the geometries. In such way, Klein presented 
9 two-dimensional spaces. However, 6 of them he considered practic unaplicable [1]. Till 
now speaking about "non-euclidean geometry", elliptic or hyperbolic geometry is primarily 
understood. Obviously, in order to make all geometries to be taken seriously, an accessible 
model is required. One of such model for two-dimensional case proposed [ I] Isaak Moiseevich 
Yaglom (1921 — 1988), using the notion of generalized complex number. Among more recent 
results you can refer to [5, 6, 7]. 

In this work, supposed to your attention an uniform model of geometric spaces and based 
on it general analytic geometry are described. Among its advantages there are its universality 
and linearity, hence easyness to use. It isn't limited to specific dimension. 

The first chapter describes different types of distance and angular measure and their mod- 
els. Different variants of axioms valid for different geometries are analyzed, as well as one 
variant of them, depending on some parameter and universally valid. A analytic model de- 
pending on some parameters is constructed. Lengths and angles are defined as parameters of 
corresponding motions. 

In the second chapter you can find triangle equations valid for all geometries. The chapter 
describes generalized orthogonal matrix as general form of motion matrix. A vector approach 
will be shown for description of points, lines and planes, and for linear calculus of lengths and 
angles. At the end of chapter, the reader will find a linear way to calculate volumes. 

The third chapter has more philosophical character then practical one. Your attention 
will be set on proper terminology and several well known spaces will be described in terms of 
constructed theory. 

Uniform model of geometric spaces becomes the background of the GeomSpace project 1 . 



1 http://sourccforgc.nct/pi'ojects/geomspace/. The last version of this book can also be downloaded form 
here. 



Chapter 1 



Geometric Space Model Construction 



1.1 Three Kinds of Plane Rotations. Rotation Charac- 
teristic 

Consider real plane M 2 . Consider three different transformations of IR 2 : rotation 9V (</>), 
Galilean transformation 9t"(</>) and Lorintz transformation 9V" (0) defined by matrices: 



cos <p — sm ( 
sin 6 cos </> 



1 



and 



tfliir ( L\_ /cosh0 sinh<A 
ysinh0 cosh</>y 



where <fi G IR. 

Transformations 9V(0), 9V'(0) and 9t"'(0) have several common properties. The determi- 
nant of all their matrices is 1, all them have the only fixed point — origin O = (0, 0), 91(0) = J 
- unit matrix, 9t(x)9t(?/) = 9t(x + y) = %\(y)*R(x) and the trajectory of point P = (1, 0) 
verifies equations: 

x 2 + x \ = 1, for 91', 

xo = 1, for 91", 
xl — x\ = 1, for 91'". 
More general, the trajectory equation can be written as (Figure 1.1): 

0Cq ~\~ li — — X ^ k — — X ; ; X . 

We will name 91' elliptic rotation, 91" parabolic rotation, 91"' hyperbolic rotation and 
respective angle. We will name the coefficient characteristic of a rotation, k = 1 corresponds 
to elliptic, k — to parabolic and = — 1 to hyperbolic rotation. 
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Figure 1.1: Trajectory of point P on transformations fK', 9t" and 9V". 



1.2 Functions C(x). Six) and Tlx 



We can see, that the matrices 9V , 9V' and 9V" have elements rn = and ri2 = —kr2i- We 
can write: 

where 

{cosx, k — 1 

1, fc = 

coshx, fc = — 1 

and 

sinx, k = I 

S(x) = { x, k = 

sinhx, k = —1 

Finally, we can define formally the functions C(x) and S(x) as 1 : 

C(x) = C(x,k) = J2hkT-—, (1.2) 

n=0 

p 2n+l 

(2n+ 1)!' 

n=0 V ' 

We will introduce one more function: 



S(x) = S(x,k) = fl(-kr-^— Vl . (1.3) 



Note, that always has place the equality: 

Here and further we will consider for simplicity that k — 1 for fc = too. We will say a: divide fc 1 , fc = 
if in expression x/fc z the exponent of fc in numerator is not less then i. 



1 .3. REPRESENTATION OF TRANSLATION AS ROTATION. ITS CHARACTERISTICS^ 



C 2 (x) + kS 2 (x) = 1, Vsc G K. (1.5) 

We will name transformation (1.1) generalized rotation. Note, that along the angle it has 
one more parameter — its characteristic. 

1.3 Representation of Translation as Rotation. Its Char- 
acteristics 

Generalized rotation has a fixed point. However, the translation usually doesn't have a fixed 
point 2 . We will define the translation through generalized rotation using an extra-dimension, 
as it is done in projective geometry. 

For n-dimensional space consider vector space R n+1 . Rotation matrices have two different 
rows compared to unit matrix. When one of them is the first row, we will consider them 
translations. When none of them is the first one, we will consider them rotations. Therefore, 
the first coordinate (we will count it 0) will be additional. 

Let n — 1. We will name the vector o = {1,0} origin. If a = 9t(</>)o, then the angle 
between o and a we can name distance oa. Different values of k correspond to different 
translations: when k = 1 translations are elliptic, when k = they are parabolic, and when 
k = —1 they are hyperbolic. 

There is a kind of distance measure for each kind of translation: elliptic, parabolic and 
hyperbolic. The difference between them can be seen in variants of V postulate of Euclid for 
elliptic, linear and hyperbolic geometry (Figure 1.2). 




a) b) c) 



Figure 1.2: Variants of Euclid's V postulate — elliptic a), linear b) and hyperbolic c). 

Elliptic postulate (Figure 1.2 a) is 3 : For a given line I and a point P ^ I, exists no line 
p 3 P 1 1 Dp = 0. It is identical to the following: For a given line / and a point P ^ /, all lines 
p 3 P intersect /. 

The linear postulate (Figure 1.2 b) is: For a given line I and a point P ^ I, exists one line 
p 3 P | / n p = 0. 

The hyperbolic postulate (Figure 1.2 c) is: For a given line I and a point P I, exist at 
least two lines j/, p" 3 P 1 1 D j/ = 0, / H p" = 0. 

2 A fixed point is present for example in translation on elliptic plane. 

3 For this case it is necessary to modify another two postulates, namely that from any three points on a line 
exactly one lies between two others, and that any line can be extended infinitely in any direction. 



10 



CHAPTER 1. GEOMETRIC SPACE MODEL CONSTRUCTION 



Generally, V postulate of Euclid can be formulated as: For a given line / and a point P ^ I, 
exist fcl lines p 3 P \l fl p = 0. It should be mentioned that fcl is a symbol, not a number 
used in calculus. Its value equals to for k\ = 1, 1 for k\ = and oo for ki = — 1. 



1.4 Kinds of Space Rotations. Bundles of Unconnectable 
Points 

It's easy to see that classic rotations in Euclidean geometry, as well as in the elliptic (Rie- 
mannian) geometry and the hyperbolic (Bolyai-Lobachevsky) geometry has the characteristic 
k — 1. We can extend the notion of space rotation to generalized space rotation with some 
characteristic. The best way to illustrate difference between them is to formulate angular 
equivalent of V Postulate of Euclid — axiom of points connectability (Figure 1.3). In order 
to do this we will change the following phrases between them: 

line I < — > point L, 

P e I < — > p3 L, 

P(£l < — ► p^L, 

AB = (j) < — > Zab = 0, 

a n b = C < — > c = AB, 

a(lb = < — > A is unconnectable with B. 



The last statement is unusual for the above three geometries . It makes sense in geometries 
with angular characteristic or —1. The unconnectable property of points is similar to parallel 
property of lines. 




Figure 1.3: Different variants of points unconnectability axiom — elliptic a), linear b) and 
hyperbolic c). 



The angle equivalent of V Postulate for elliptic characteristic (Figure 1.3 a) is: On a line 
/ ^ P exist no points L unconnectable with P. 

For parabolic characteristic (Figure 1.3 b) it is: On a line I ^ P exists the only point L 
unconnectable with P. 

For the hyperbolic characteristic (Figure 1.3) it is: On a line / ^ P exist at least two points 
V and L" unconnectable with P. 



It conflicts with axiom which states that through any two points goes a line. This axiom should be changed 
by one of the following in order to consider the geometries with non-elliptic rotations. 



1.4. KINDS OF SPACE ROTATIONS. BUNDLES OF UNCONNECTABLE POINTS 11 



Generally this axiom can be formulated as: On a line I P exist k2 points L unconnectable 
with P. As in case of parallel lines, symbol k2 isn't used in calculus. 

Similar to bundles of lines — intersected, parallel or divergent we can speak about bundles 
of points. More exactly, let X,Y £ lR n+1 . All linear combinations Z = aX + /3Y, a, (3 £ R form 
a set we will name bundle of points. As we will see, this set has one constraint. Therefore, 
it has one free parameter. As every two lines define a bundle of lines, every two points (X 
and Y) define bundle of points. If X is connectable with Y this bundle is a line (similar to 
intersection point of bundle of intersected lines). Lines has blue color on figure 1.3. If X 
and Y are unconnectable, this bundle of points isn't a line (similar to bundle of parallel or 
divergent lines). Bundles of unconnectable points are green and red on figure 1.3. 

For any angle characteristic there are infinity of bundles of connectable points. For angle 
characteristic 1 all point bundles are lines. For angle characteristic for any point there is 
the only bundle of unconnectable points (green). For angle characteristic —1 there are infinity 
bundles of unconnectable points (red). In thes case the bundles of connectable points and 
the bundles of unconnectable points for some point form two categories of bundles. The limit 
(marginal) bundles of unconnectable points (green) can be viewed as the third category (similar 
to differencee between parallel and divergent lines). There are exactly two limit bundles. Note 
that bundles of connectable points intersect all circles with centre in the centre of bundle, all 
bundles of unconnectable points don't intersect these circles and limit bundles are asymptotic 
to circles (Figure 1.4). 




a) b) c) 

Figure 1.4: Mutual position of different bundles and circles a) elliptic angular characteristic, 
b) linear angular characteristic and c) hyperbolic angular characteristic. 



Emphasize that the angle between two lines and the angle between two two-dimensional 
planes are the different measures. The angle between two threedimensional planes is different 
from them both and so on. Thus, the angle between lines can have the different characteristic 
then the angle between two-dimensional planes and so on. 
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1.5 Main Space Rotations 

Consider R n+1 and k u k 2 ,...k n G {-1,0,1}. We will note d(x) = C(x,ki), Si{x) = S{x,ki) 
and Ti(x) = Si(x)/Ci(x). Let 
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We will name fHi, ...!SK n mam space rotations. 



1.6 Vector Product. Invariant Quadric Form 



Let 



K m = Yl h, Vm = 0, 



n 



i=l 



(1.6) 



We can see that X m G {—1,0, 1}, Vm = 0, n as well as k m . Let define vector product as 

n 

xQy = s ^K i x i y i (1.7) 

i=0 

For some vectors x = {x , x ± , ...x n } and y = {y , y ± , ...y n }, x' = %n(<f>)x = {x , ...x m _ 2 , x m _i 



1 . 7. SPACE DEFINITION BY ITS SPECIFICATION 
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S m {(p) + x m C m (<p),x m+1 , ...x n } and y' = ^K m {4>)y. We can see that 

n 

x'Qy' = ^Kix'-y't 

m-2 

= ^ K i X iUi 
i=0 

+ \\ x m-\Cm\4>) ~ k m X m S m {(f>)){y m -iC m [(p) — k m y m S m ((fi)) 

+ k m\Xm— 1 S m (<f>) + x m C m ((j)))(y m -iS m ((j)) + y m C m ((p)))K m -i 

n 

+ KiXiyi 

i=m+l 
m-2 

i=0 

~l~ ((^m— lVm— l^m(^) k m {x m —\y m -\- X m y m —\ )S m {4>)C m (4>) + k^XmymS^^) 
m%m—iym—iS m ((fi) + k m {x m —\y m 4~ x m y m —\ 

)Sm{4>)Cm{4>) + k m X m y m Clj^>y))^m-\ 

n 

j=m+l 
m-2 

i=0 

+ (x m _iy m _i(C^(0) + /c m S^(0)) + /c m a;ml/m(C^(0) + /E m S^(0)))if m _i 

n 

i=m+l 



^ = xQy 



i=0 



This is true for all m = l,n. So the quadric form 10?/ is invariant in respect to main rotations 

Of Vtn. 



1.7 Space Definition by its Specification 

Consider MP" projective space and fc, G {— 1,0, 1}, Vi = l,n. We can now introduce a 
geometric space 'unit sphere' B" = {1 6 MF n \ x Q x = 1} (Figure 1.5). As all main rotations 
preserves the quadric form defined by product 0, they also preserves B n . We will name 
ki, i — 1, n space specification. We will name 'point' XeB" the corresponding vector x G lRP n 
and will use homogeneous coordinates normalized in order to x x — 1. 

We will name 'origin' of B n the point O = [1 : : ... : 0] G l n . It isn't origin of R n+1 , 
(0, 0, ...0) ^ B n and we will refer to O as origin if isn't specified otherwise. 

It's easy to see that for any k x , fc 2 , ...k n , O = 1° C I 1 C ... C B n . 

We will define motions of B n all transformations that result on finite product of main 
rotations. 
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Figure 1.5: Sphere of space with specification { — 1, — 1}. 



We will define 'lines' all images of B 1 on any motion of B n . Similarly, we define 'Tri- 
dimensional' planes all images of B m on any motions of B n for any m G 0, n — 1. 

For each characteristic parameter ki we can introduce a scale parameter e K+, i = l,n. 
The k\jr\ is exactly the gaussian curvature of space. Others have no representation since 
finite angle measure doesn't require scaling. In this case the radian measure is native. An 
example of angle scale is degree measure which has scale 180/7T. However when the angle is 
not bounded a scale introduction has sense. All scales can be easy embedded in functions 

Ci(x), Si(x) and T t (x) by using instead d m, $ (fA and T { (fA respectively, i = l~h. 



1.8 Definition of Measure Using Motions 

A traditional way of definition the measures and motions is to provide a way to calculate the 
distances as is and then to define motions in such way that all maps M : W 1 — > W 1 preserve 
the distance. We go another way. We provide motions as is and then search for a way to 
define measures in such way that motions preserve them. 

We will say point A e B 1 C B n has the distance from origin O if A = *Ri(<j))0. Having 
O = [1 : : ... : 0], A = [d(0) : : : ... : 0], 0&A = We will say one-dimensional 

(planar) angle between B 1 and some one-dimensional line B' 1 C B 2 equals <fi if B' 1 = ^(^B 1 . 
Similarly, we will define the m-dimensional angle <fi between B m and m-dimensional plane 
B /m C B m+1 if M' m = fH m+1 (0)B m , Vm = 0,m - 1. Note, that n-dimensional angle between 
any planes is since all them are subset of B n . 



1.8. DEFINITION OF MEASURE USING MOTIONS 
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Let X, Y e B n . If there exists a motion that maps B 1 to XY we will name points X and Y 
connectable and distance XY measurable. If not, we will name points X and Y unconnectable 
(just as lines can be parallel) and strictly speaking the distance XY doesn't exists 5 . 

We can find a motion 501 of space M n that maps origin O to X and some point A e B 1 C B n 
to Y. As motion 9Jt preserves the quadric form 0, we can see that X Y = O A. We can 
define the distance <fi between X and Y as 



It's easy to see that all motions preserve the distance. In case of elliptic, Euclidian and 
hyperbolic space it is sufficient, because all other measures can be calculated from distances. 
However, in some spaces angles can be scaled in a manner distances are scaled in Euclidean 
space. So we should find the way to measure all the measures in general case. 



5 In this case there exist a measure XY, but it may have different characteristic then distance. We will 
name this measure also distance, keeping in mind that it is generalized distance. 



d(0) =XQY. 
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Chapter 2 
Measure Calculus 



2.1 General Triangle Equations 

Consider triangle ABC £ B 2 with the edges a, b, c, interior angles a, 7 and exterior angle 
P' (Figure 2.1). Let A = [1:0:0] the origin, C = 5Ki(6)A = [Ci(6) : Si (6) : 0] and 5 = 
9 c l2(a)9 c ii(c)v4 = [Ci(c) : 5i(c)C2(a) : 5i(c)S f 2(a)]. Note, that the interior angle (3 does not 
exist in case of k 2 = or k 2 — — 1. The exterior angle /?' always exists. 




Figure 2.1: General triangle equations deduction. 



Now, let A'B'C = 9l(—b)(ABC) (Figure 2.1, cyan). The point we are interested in is 
B> = [Ci(6)C 1 (c) + k x Sx{b)S x {c)C 2 {a) : -Si(6)Ci(c) + C 1 (6)5 1 (c)C 2 (a) : Si(c)S 2 (a!)]. At the 
other hand, now B' = [C\{a) : — Si(a)C2( 7 ) : Si (0)5*2(7)]. It means: 

d(o) = C 1 {b)C l {c)+k l S 1 {b)S l {c)G 2 {a), 
-5 1 (a)C 2 ( 7 ) = -S^C^c) + C 1 (6)5 1 (c)C 2 (a), 
Si(a)S 2 ( 7 ) = Si(c)S 2 (a). 

The first equation is the form of the Cosine I law. Similarly we have 

d(c) = CMCiib) + fei5 , i(o)5 , 1 (6)C 2 ( 7 ). 
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CHAPTER 2. MEASURE CALCULUS 



The third equation is equivalent to 

51(a) Stic) 



S 2 (a) 52(7)' 

which is the form of the Sine law. 

Let now A"B"C" = 9ti(-c)0ta(-a)(AB(7) (Figure 2.1, brown). Now we are interested in 
vertex C" = [Ci(6)C7 1 (c) + fei5'i(6)5 , i(c)C7 2 (a) : -Ci(&)o\(c) + o\(&)Ci(c)C 2 (a) : -o\(fc)S 2 (a)]. 
At the other hand, C" = [C?i(a) : o\(a)C 2 (/3') : -S 1 (a)5 2 (/3')]- From here we have: 

d(a) = C 1 (6)C 1 (c) + A;iS , 1 (6)S , 1 (c)C 2 (a), 
5i(a)C 2 (/3') = -Ci(6)5 , i(c) + 5'i(6)Ci(c)C 2 (a), 
-S^a)^) = -S 1 (b)S 2 (a). 



The first one is the Cosine I law, the third one is equivalent to: 

ffi(q) = Sijb) = gi(c) 
52(a) W) S 2 ( 7 ) 



(2.1) 



which is the Sine law. Note that in case k 2 = 1 we have /3 = 7r — /?', S 2 (f3) = 5 2 (/3'). Let 
calculate the value of C 2 (a) from the first equation and put it to the second one: 

sMcm = -c^ic) + Stib^icf 1 ^ ~ Clib)Cl(c) 



-C^b) Stic) + 0,(0) 



hS^Stic) 
Ci(a)-Ci(6)Ci(c) 



hSi(c) 

hS^S^C^') = -hS^cfCt^ + C^C^-Ctib^c) 2 
= Ci(a)Ci(c) - d(6)(Ci(c) 2 + fciSi(c) 2 ) 
= CtiajC^-Ctib), 
d(6) = ^(^^(^-^^(^^(c)^^'). 

Note the ' — ' sign in the right part of the equation. It is so because the (3' angle is external. 
For the case k 2 = 1, the internal angle (3 = ir — (3', C 2 (/3) = — C 2 {j3'). 
What about the Cosine II law? We will use these two equations: 

-o\(a)C 2 ( 7 ) = -5' 1 (6)C 1 ( C ) + C 1 (6)5'i(c)C2(q ; ), 
5i(a)C 2 (/3') = -Ci(6)5 1 (c) + 5i(6)C7i(c)C 2 (a) 

First, replace Si (6) with Si(a)S 2 (/3')/S 2 (a) and Si(c) with Si(a)S 2 ( 7 )/S 2 (a): 

-5i(g)C 2 (7) = -5i(a)^C' 1 ( C ) + C'x(6)5 1 (o)f^C' a (a), 

o 2 (aj 02(a) 

-S 2 (a)C 2 ( 7 ) = -C 1 (c)S , 2 (/3 , ) + ^i(&)5 2 (7)C 2 (a), 
S2(P')C 1 (c) = S 2 (a)<7 2 ( 7 ) + C 2 (a)S 2 ( 7 )Ci(&), 



2.1. GENERAL TRIANGLE EQUATIONS 
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and 

Sx{a)C 2 tf) = -C 1 (b)S 1 (a)^ + S 1 (a)^flc 1 (c)C 2 (a), 

02(a) 02(a) 

S 2 (a)C 2 ((3') = -C 1 (6)5 2 ( 7 ) +C 1 {c)S 2 {F)C 2 {a), 

5 2 (7)Ci(6) = -5 2 (a)C 2 (/3') + C 2 (a)5 2 (/3 , )Ci(c). 

Now from the first equation let calculate Ci(c) and put it in the second one: 

<? (~\n (h\ q ( m\c ( r'\ 1 n ( n \Q (ot\ S?(PdCM + C 2 (a)S 2 (-f)d(b) 
^2(7)^1(0) = -02(a)C 2 (p ) + C 2 (a)o 2 (p ) S^jf) 

= -S 2 {a)C 2 {p) + C2(a)S2(a)C 2 (i) + C 2 (a) 2 5 2 ( 7 )C 1 (6), 

1 S 2 ( 7 )C 1 (6)(l-C 2 (a) 2 ) = S- 2 (a)(C 2 (a)C 2 ( 7 )-C 2 (/3')), 

fc25 2 (7)Ci(6)5 2 (a) 2 = S 2 (a)(C 2 (a)C 2 ( 7 )-C 2 O0)), 

fc25 , 2 (a)S , 2 ( 7 )Ci(6) = C 2 (a)C 2 (7)-C 2 (/3'), 

W) = C 2 (a)C 2 ( 7 )-A; 2 5 2 (a)5 2 (7)C 1 (6). 

When fc 2 = 1 we have 

— cos (3 = cos a cos 7 — k 2 sin a sin 7 Ci(fe), 
cos/3 = — cos a cos 7 + k 2 sin a sin 7 Ci(6). 

Similarly, calculating C\{b) form the second equation and putting it in the first one, obtain: 

S 2 {P)C 1 (c) = S 2 (a)C 2 ( 7 ) + C 2 (a)S 2 ( 7 ) — 

= S2(a)C 2 (i) + C 2 (a) 2 S 2 (0')C 1 (c) - C 2 (a)S 2 {a)C 2 {p), 

5 2 (/3')Ci(c)(l-C 2 (a) 2 ) = 5 2 (a)(C 2 ( 7 )-C 2 (a)C 2 (/3')), 

foSitf'McWa) 2 = 5 2 (a)(C 2 ( 7 )-C 2 (a)C 2 (^)), 

k 2 S 2 (a)S2^)C 1 (c) = C 2 ( 7 )-C 2 (a)C 2 (/3'), 

C 2 ( 7 ) = C 2 (q ; )C 2 ( / 5 / ) + A; 2( 5 2 (q ; ) ( 5 2 ( / 5 / )C 1 ( C ). 

When fc 2 = 1 we have as above 

cos 7 = — cos a cos /3 + A; 2 sin a sin /3Ci (c) . 

Similarly, we have 

C 2 (a) = C 2 (/3')C 2 ( 7 ) + fc 2 S 2 03')5 2 (7)Ci(a). 

We will find the form of the Cosine I and II law that does not contain C\ or C 2 functions 
in the left part. However, it contains these functions in the right part. It makes sense since 
in the case k\ 7^ (for the Cosine I law) and k 2 7^ (for the Cosine II law) when we can 
calculate their respective C^ 1 functions, but when k\ = 0, the space admit distance scaling and 
the angle values does not determine the distances (Cosine II law is a equality which doesn't 
contain C\ function), while when k 2 = 0, the space admit the angular scaling and distances 
does not determine angles (Cosine I law is a equality which doesn't contain C 2 function). 
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Note also that we can deduce one form of Cosine I and one form of Cosine II law if we 
introduce a (may be virtual) angle j3 so as: 



S 2 (p) = s 2 (J?), 

cm = -W), 

T 2 (P) = -UP 1 ). 



Then both Cosine I and II law have identical form. Now let calculate 



fciS?(a) = l-Cf(a) 

= (01(b) + hSfmClic) + hSl(c)) 

- (C7i(6)Ci(c) + fciS'i(6)S , i(c)C 2 (a)) 2 

= Cl(b)Cl(c) + fc!(^(6)^(c) + k^(b)dl(c) + fc 2 S 2 (&)S 2 ( C ) 

C{(b)Cl(c) - 2fciCi(6)Ci(c)S , i(6)S , i(c)C 2 (a) - k*S*(b)S*(c)C*(a) 

= fc^C^Sftc) + Sl(b)Cl(c) - 2C 1 (b)C 1 (c)S 1 (b)S 1 (c)C 2 (a)) 

+ fc 2 S 2 (&)S 2 (c)(l-C 2 2 (a)) 

= fc 1 ((^(6)S' 1 2 (c) + Sl(b)C*(c) - 2C 1 (b)C 1 (c)S 1 (b)S 1 (c)C 2 (a)) 

+ fc?fc 2 S?(&)S?(c)Si(a), 

S?(a) = C 1 2 (6)5 1 2 (c) + 5 1 2 (6)C 1 2 (c)-2C 1 (6)C 1 (c)5 1 (6)5 1 (c)C 2 (a) 

+ fc 1 A; 2 ^(6)S' 1 2 (c)Sf(a), 



or, having: 



Similarly, 



and 



d(a) = Ci(6)Ci(c)(l + fciri(6)T 1 (c)C 2 (a)), 

T 2 ( v gg) + Tflc) - 2T 1 (6)r 1 (c)C 2 (a) + hhTfWTftfSfta) 
11 J (l + fc 1 T 1 (6)r 1 (c)C 2 (a)) 2 • 1 • J 



5 2 (6) = C 1 2 (a)5 1 2 (c) + 5 1 2 (a)C 1 2 (c) + 2C 1 (a)C 1 ( C )5 1 (a)5 1 ( C )C 2 (/3 / ) 

+ A; 1 A; 2 5 1 2 (a)5 1 2 (c) 1 S 2 (/3 , ), 

T 2, w T 2 (a) + T 2 (c) + 2T 1 (a)T 1 (c)C 2 (/3 / ) + hhT^T^SfiP') 
l{ ) " (1 - fciT^T^CW')) 2 ~ " ' ( j 



5 2 (c) = C 1 2 (a)5 1 2 (6) + 5 1 2 (a)C 1 2 (6)-2C 1 (a)C 1 (6) 1 S 1 (a) 1 S 1 (6)C 2 (7) 

+ h 1 k 2 %(a)Sl(b)S*(>y), 

rytii \ Tf(a) + Tl(b) - 2T 1 (a)r 1 (fr)C 2 ( 7 ) + fc 1 fc 2 r i 2 (a)r i 2 (6)5 1 2 ( 7 ) 
lU " (l + fc 1 T 1 (a)T 1 (6)C 2 ( 7 )) 2 ' 1 ' ] 
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Now, let calculate 

k 2 S 2 2 (a) = l-C 2 2 (a) 

= (C 2 2 ((3') + k 2 S 2 2 ((3'))(C 2 2 (l) + fcSf (7)) 

- (C 2 (^)C 2 ( 7 ) + hS^S^C^a)) 2 

= C 2 2 (/3')C 2 2 (7) + hCK^S 2 ^) + k 2 S 2 2 (P')Ci (7) + fc 2 2 5 2 2 (/3')5 2 2 ( 7 ) 

C 2 (/3')^ 2 2 (7) - 2fc 2 C 2 (/3')5 2 (/3 , )C2(7)^(7)Ci(a) - *|Sf (/^(^(a) 

= A; 2 (C 2 2 (/3')S 2 2 ( 7 ) + 5 2 (/3')C|(7) - 2C 2 {p')S 2 {P')C 2 { 1 )S 2 { 1 )C l {a)) 

+ k 2 2 S 2 2 ^)S 2 2 mi-Cl{a)) 

= k 2 (C 2 (P')S 2 M + S 2 (P')C 2 M - 2C 2 ((3')S 2 ((3')C 2 (~f)S 2 (~f)C 1 (a)) 

+ M 2 S 2 2 (/30S 2 2 ( 7 )S 2 (a), 

S 2 2 (a) = C 2 2 (/3 / )5 2 2 (7) + 5 2 (/3 / )C 2 2 (7)-2C 2 (/3 / )5 2 (/3 / )^ 2 (7)52(7)Ci(a) 



or 

T 2,„ x _ T|(/3Q + T|( 7 ) - 2r 2 (/3Qr 2 (7)C 1 (a) + k^T^T^Slja) 



Similarly, 



S 2 (P') = C 2 («)5 2 (7) + S 2 (a)C 2 (7) + 2C 2 (a)^ 2 («)C 2 (7)^(7)^iW 

+ hhsl^sl^slib), 

r2(R n _ T 2 (a) + T|( 7 ) + 2^(0)^(7)^(6) + fc 1 fc 2 r 2 2 (a)r 2 2 (7)5 1 2 (6) 
aVU " (l-A^a)^)^)) 2 ' {I - b) 

and 

S 2 2 ( 7 ) = C 2 («)5 2 (/30 + 5 2 («)C 2 (/30-2C' 2 (a)S 2 (a)C 2 (/3 / )5 2 (/3 , )^i(c) 
+ hhSii^Sii^SKc), 
2 = T 2 (g) + T 2 (f3') - 2T 2 (g)T 2 (/3QC 1 (c) + fcifc 2 T 2 2 (g)7f (^ ; )5?(c) 
2l7J (l + A; 2 T 2 (a)T 2 (/3')Ci(c)) 2 " 1 ' J 

What does it mean for triangle? From the Sine law (2.1), having function S(x) monotoni- 
cally increasing result that the longest side of any triangle is opposite to the largest angle and 
the shortest side is opposed to the smallest angle. From the Cosine I law in its form that uses 
C\(x) function, having 

C a (6)Ci(c) + hS^Stic) = d{b- c), 

Ci(a)Ci(6) - fci^i (0)51(6) = Ci(a + 6), 

and Cj(a;)Zel and is decreasing when ki = 1, Cj(x) = 1 and is constant when ki = 0, Ci(x)gel 
is increasing when ki = —1, we can see: 

d(a) = Ci(&)tfi(c) + fei5 1 (6)5 1 (c)C 2 (a) 
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is equivalent to 





::: 


c, 


A-2 


= 1 


a | 




c, 




= 




[<&- 


c, 


h 


= -1 



Similarly, 

d(6) = C 1 (a)C 1 (c) - fe a S' 1 (a)5 1 (c)C a (/3') 

is equivalent to 

{< a + c, £; 2 = 1 
= a + c, fc 2 = 
> a + c, &2 = — 1 

From the Cosine II law we can see: 

C 2 (a) = C 2 (^)C 2 ( 7 ) + k 2 S 2 {f3')S 2 { 1 )C 1 {a) 

is equivalent to 

r > p - 7 , h = i 

a J = (3' - 7 , fci = 
7 , fc 1 = -l 

Similarly, 

C 2 (^) = C 2 (a)C 2 ( 7 ) - A; 2 5 2 (a)5 2 ( 7 )C 1 (6) 

is equivalent to 

< a + 7, fci = 1 
>" { = a + 1: h = 
> a + 7, fci = — 1 

2.2 Right (Quasi)— Triangle Equations 

We can define orthogonality in B n using the orthogonality in IRP n . Namely, two vectors V\ 
and v 2 of space 1RP™ are orthogonal, if v±Q v 2 = 0. 

For B 2 plane and a line, the orthogonal bundle is line only if k 2 — 1. In this case when line 
rotates count-clockwise, its orthogonal line rotates count-clockwise and vice-versa (Figure 1.4 
a). When k 2 = there is the only orthogonal bundle, which doesn't rotate (Figure 1.4 b). 
When k 2 = — 1 the orthogonal bundle rotates clockwise when the line rotates count-clockwise 
toward to the same limit bundle and vice- versa (Figure 1.4 c). 

Generally we can't speak about right triangle as one of its catheti is line and another 
isn't (when k 2 ^ 1). However, as we will see, this figure is important. We will name it 
right (quasi)-triangle, which means right triangle, when k 2 = 1 and right quasi-triangle when 
k 2 / 1. 

We will construct a (quasi)-triangle as half of isosceles one (Figure 2.2). Consider a triangle 
A B A' with A Bq = A' Bq = c, A A' = 26, ZA' A B = ZA A' B = a and external angle 
ZA B A' = 2{3'. 



2.2. RIGHT (QUASI)-TRIANGLE EQUATIONS 

B B Q 
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Figure 2.2: Right (quasi)-triangle equations deduction. 



Let A' = O = [1 : : 0] be origin, A = VKx{2b)A' Q = [Ci(26) : Si (26) : 0], B = «K 2 (a)9ti(c) 
A = [Ci(c) : S x {c)C 2 {a) : S 1 (c)S 2 (a)}. 

Let ABA' = ^ 1 (-6)(A 0j B A' ) (Figure 2.2, black). Now, A' = 9li(-&Ko = l C i( b ) : ~ S i( b ) 
0], A = 9ti (-6) A, = [Ci(6) : 51(6) : 0] and B = JKi(-&)fl = 



d(6) hSiQ)) N 
-^(6) d(6) 
1. 



C x (c) 
Si(c)C 2 (ar) 



= [Ci(6)Ci(c) + fc 1 S' 1 (6)S , i(c)C 2 (a) : -Si(6)d(c) + d (6)5i(c)C 2 (a) : S 1 (c) S 2 (a)}. 

Finally, let C E AA', AC = A'C = b. Then C = [1 : : 0] is origin. From figure equality 
A'BC = ABC result BC _L A' A. Therefore we can consider ABC right (quasi)-triangle. 

Having figures A'BC = ABC and C is origin, result B have form B = (x, 0, y), where 

dtydic) + hS^S^Czia) = x, 
-S 1 (b)C 1 (c) + C 1 (b)S 1 (c)C 2 (a) = 0, 
S 1 (c)S 2 {a) = y. 

From the second equality, have 

T 1 (b) = T 1 (c)C 2 (a). (2.8) 
Using the value of C 2 (a) from this equality and putting it in the first one, have 

Ti(b) 



x 



C 1 (b)C 1 (c)+k 1 S 1 (b)S l (c) 
CAc) 



d(6) 



(Cf^+hSiib)) 



Ti(c) 

Ci(c) 

Ci(6)" 
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Let now calculate the value of 

x 2 + kik 2 y 2 = + kik 2 Sl(c)S%(a) 



cm 

cm 



^ + k 1 S 2 (c)(l-C 2 (a)) 

|| + ^( C )-^( C )|| 
Cf{c) + hS 2 (c) - hSl{b) Cl{c) 



C 2 (b) 1 1W 1 1V ' C 2 (b) 
-(1 - kiSf{b)) + k x S 2 {c) 



C 2 (b) 
C!(c 



C 2 (b) + hSUc) 



Ci(c)T 2 (a), 



C 2 (6) 

= C?(c) + A;iS?(c) = l. 

It means that exists a G K, Ci 2 (a) = ^r^y,Si 2 (a) = Si (0)62 (a) that has characteristic A; = 
k\k 2 . It is a 'distance' parameter BC . We have two more equations: 

d(c) = C 12 (a)d(6) (2.9) 
Si 2 (a) = 51(0)5-2(0) (2.10) 

From (2.8) and (2.10) have: 

5i 2 (a) _ S 1 (c)S 2 (a) 
Ti(6) ~Ti(c)C 2 (a) 

using the value of Ci(c) from (2.9) have 

^ = Ci2(G)Ci(6)T 2 (a), 

T 12 {a) = S 1 (b)T 2 (a). (2.11) 

The last 6 equations will include (3'. In order to be able to deduce them we will introduce 
translation T(— a), so as T(— a)B = C. Having characteristic a is k\k 2 = K 2 , 

( C 12 (a) K 2 S 12 (a) 
%{-a) = 1 

\-S 12 (a) Ci 2 (a) 

We can check this map preserves vector product. 

Applying T(-a), obtain B l = %(-a)B = [1 : : 0], C x = %{-a)C = [C 12 (a) : : -5 i2 (a)] 
and A x = %(-a)A = 

C 12 (a) K 2 S 12 {a) 

1 
S 12 (a) C 12 (a) 
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= [C 12 (a)C7 1 (6) : S 1 (b) : - S 12 {a)C 1 (b)} = [C x {c) : S x {c)C^') : -^(c)^/?')] (Figure 2.2, cyan). 
From here we have 

S 1 (b) = S 1 (c)C 2 {p). (2.12) 

Moreover, having (2.9), obtain 

S 12 {a)C 1 {b) = S 1 (c)S 2 ^'), 

SM^Q = s 1 (c)s 2 (n 

T 12 (a) = T^S^P'). (2.13) 

Combining the last 2 equalities (2.12), (2.13) with (2.9), we have 

T 12 (a) Ti{c)S 2 {P) 

Si(b) S x {c)C 2 (P') 

W) T 2 (/3') 

~ Ci(c) ^ Ci^C^)' 

S 12 (a) = T x {b)T 2 {P'). (2.14) 



Now, having (2.8) and (2.12): 

T 1 (c)C 2 (a)=T 1 (6) 

calculate with (2.9): 

Ci(6) = 



gig) _ gi(c)C 2 (^) 
" Ci(&) ~ Ci(&) 

gi(c)C 2 (/3Q 
T 1 (c)C 2 (a) 



C 12 (a)d(6) 



C 2 (a) 

C 2 (/?') 



C 2 fal 



C 2 (a) = C 12 (a)C 2 (/3'). (2.15) 

Now from (2.11), (2.12) and (2.13), 

T 12 (a) = 5 1 (6)T 2 (a) = T^c)^ (/?'), 
S 1 (c)C 2 (/3 / )T 2 («) = Ti(c)S 2 {p), 

T 2 (/3') = d(c)T 2 (a). (2.16) 

Finally, by multiplying the last equations (2.15) and (2.16), have using (2.9): 

C x {c)T 2 {a)C 2 {a) = T 2 {f3')C 12 {a)C 2 {/3') , 

Ci(c)5 2 (a) = C 12 (a)C 2 (l3% 

Ci 2 (a)C7i(6)5 , 2 (a) = C 12 (a)S 2 (f3'), 

S 2 (f3') = d(&)S 2 (a)- (2.17) 
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It is necessary to modify equations (2.9) and (2.15) in order to not contain the C(x) 
function. 

fciSf (c) = 1 - C 2 (c) = (C 2 2 (a) + hk 2 S 2 2 (a))(C 2 (b) + hS 2 (b)) - C 2 l2 {a)Cl{b) 
= hk 2 S 2 2 (a)C 2 (b) + hCUa)S 2 (b) + klk 2 S 2 12 {a)S 2 x {b), 
S 2 (c) = k 2 S 2 12 (a)C 2 (b) + C 2 2 (a)S 2 (b) + k^Sf^Sfib) 

By dividing the last equality by its C(x) form, obtain: 

T 2 (c) = k 2 T 2 2 (a) + T 2 (b) + fc^l* (a)7?(6). (2.18) 

Similarly, 

k 2 S 2 2 (a) = 1 - C 2 2 (a) = (C 2 12 (a) + M 2 S 2 2 (a))(C< 2 2 (/3') + k 2 S 2 2 ^')) - C 2 2 {a)C 2 {0') 
= k 2 C 2 2 (a)S 2 2 ^') + hhSU^CKP') + hk 2 2 S 2 2 (a)S 2 2 (0'), 
S 2 2 (a) = C 1 2 2 (a)5 2 2 (/3 / ) + A; 1 5 1 2 2 (a)C 2 2 (/3 / ) + A; 1 fc 2 5 1 2 2 (a)5 2 2 (/3 / ) 

By dividing the last equality by its C(x) form, obtain: 

T 2 (a) = hT 2 2 (a) + T 2 2 (/3') + k x k 2 T 2 12 {a)T 2 2 (0') (2.19) 

Note that for k 2 = 1 equations (2.8) — (2.19) can be used if external angle 0' change to 
internal with the following changes: 



= 




cos = 


sin 0' 


sin /3 = 


cos 0' 


tan (3 = 


cot /3' 


cot = 


tan/3' 



2.3 More rotations 



As we can see, transformation %(— a) preserves vector product. In order to be a motion 
it needs to be presented as finite product of main rotations. If a, b, c, a and /3' are real 
numbers for which have place equalities (2.8) — (2.17) then it can be checked that X(a) = 
9 c l 2 (/3')9^i(c)9^2(— oe)fRi(— b). We will introduce new transformations as following: 

/l ... ... ... 0\ 






Vo 



Ci+i,...j{4>) 



Si+i,...j(<f>) 



o 



Ci+i,...j{4>) 



o 



o 



1/ 



It's easy to see that all !SKy(</>) are motions. All they can be presented as finit product of 
main rotations. We will name them rotations of the space B n . In special case, 9^(0) we will 
name them translations T,(0) of the space. 
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2.4 Generalized Orthogonal Matrix 



For WP n with given specification k p , p = l,n we will name the vector x E WP n upper i- 
normalized, i E 0, n if -^-x x — 1. For < i < j < n we will name two vectors x and y 
upper ij-orthogonal if -^-^ — x y = 0. We will name the matrix M( n+1 ) x(n+1 ) composed of 

**-rnin(i,j) 

columns q upper orthogonal if all columns q are upper i-normalized and any two columns q 
and Cj are upper ij-orthogonal. 

It's easy to see that all main rotation matrixes are upper orthogonal. Moreover product of 
two upper orthogonal matrix is upper orthogonal. Really, let X, Y are two upper orthogonal 
matrices. It means that X is composed of (xo,...,x n ) columns and Y - - from (y ,...,y n ) 
columns, and -j^-^ — Xi Xj = -^-^ — yi y 7 - = 5^ for all i,j = 0, n, where 5^ — 1, i — j and 

5ij — 0, i 7^ j. Let Z = XY with elements Zij = J2p=o x ipUpj- Let Zi and Zj be 2 columns of Z. 
Let calculate 



— ^ ] KpZpiZpj 

l*-min(i,j) ^mm(t,j) 



p=0 



1 ( if 

~T? ( X P m \Vmxi J I 

min ( l J') p=0 \mi=0 / \m 2 =0 



ri 

•Epm.2 Um.23 



-j n n n 

^ ^ ] -^p ^ ] ^ ] x pm 1 x pm2ym 1 iym2j 

min(i,j) p=Q mi= o,Tt2=0 

j n n n 

T7 ^ ^ ^ ^ ym\iVm2j ^ ^ KpXp mi Xp m2 

*min(i,j) mi=0m2=0 p=0 

n n 

min(mi,ra2)"mim2 



j n n 

^ ^ ^ ^ ymiiym2j-K m in(mi,m2)$n 



K m in(i,j) mi=0m2=0 

1 n 

"77 ^ ^ ymiym'jKn 

K min(i,j) m=Q 

$ij 



We will name the vector x E WP n lower i-normalized, i e 0,n if -^1X^=0 ~ 1- ^or 
< « < j < n we will name two vectors x and y /ower ij-orthogonal if K max ^j) Y^ p =o = 0- 
We will name the matrix M( n+ i) X ( n+ i) composed of rows lower orthogonal if all rows are 
lower i-normalized and any two rows and rj are lower ij-orthogonal. 

It's easy to see that all main rotation matrixes are also lower orthogonal. Moreover product 
of two lower orthogonal matrices is lower orthogonal. Really, let X, Y are two lower orthogonal 
matrices. It means that X is composed of (xo,...,x n ) rows and Y — from (y ,...,y n ) rows, 
where K max(iJ) £; =0 x -^f = K max(hj) ^ = for all i,j = 0^. Let Z = XY. Let 
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calculate 



n n - / n \ / n 

p=0 p p=0 p \mi=0 / \m 2 =0 



,u rt n 

Kmax(i,j) ^ ^ ^ ^ ^ ^ 2-imi X j m2 y mi py m2 p 



p=0 F mi=0rri2=0 



n n n \ 

K max {ij^ ^ ^ ^ ^ %imi% jm2 y j ~J^~ymipym2P 
mi=0ra2=0 p=0 ^ 



n n ^ 

"r?ii?Tt2 



K ma x(i,j) y j y j Xim^j!,, ,- 

mi =0m 2 =0 Amra ( rai ' m2 ) 



K ma x(i,j) ^] 



•Eim%jm 



K 

m=0 



For some upper orthogonal matrix X has place the equality 

■? i=0 
n 

E^4 = ^ 



i=0 



Let divide it to K q , q < j: 



i=0 p=l p=l 



i=0 A « i=q A « A « 



As Kg divides Kj and K i} i > q, but K q doesn't divide K iy i < q, result that for < i < j < n, 
Xij divide K q /Ki for all z < q < j, or divide Kj/Ki. 

Having for some upper orthogonal matrix X, elements Xij divide Kj/Ki, construct the 
matrix Y of the same size with elements = \J~§^Xij- The matrix Y is orthogonal one 1 (if 
may be complex, in this case it isn't unitar, but orthogonal). Really, for some i = 0,n, 

n n ^ ri 

i=0 i=0 3 3 i=0 



1 consider VO = 0. 
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and for some j\ 7^ h = 0, n, 

n 



i=0 



because x^x^ divides — ^ 32 . 



E 



o V K h K 32 



X ijl Xi 32 



K, 



K m jn(j 1 ,j 2 ) Q _ g 



X 



V)2 



For some orthogonal matrix always has place also the following equalities for %\ ^ i 2 = 0, n: 



3=0 



K 

E|4 



i=o 



3=0 



123 



E 



'A-,,A-, 2 



Xi x jXi 2 j 

3=0 j 

^mm(ii,t 2 ) / j£ X hj X i23 



K, 



max{ii,i2) 



3=0 



K-max(ii,i2) 



-0 = 



So, x^jX^j divides 



if, 



= . It means that X is also lower orthogonal matrix. 



Inverse orthogonal matrix Y~ x is easy constructed as y'^ = yij. Then 



Kj x' 
Ki ji 





Kj Xij: 



x 
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R x ij- 

Iv 3 



The last equality isn't applicable if some characteristic ki = 0. Although is true, it isn't 
determinable having the form of 0/0. If some characteristic k m = 0, m < n, the matrix has 
the form: 

'A N 
B C, 



M 



Really, for the first m columns the upper orthogonality condition is equivalent to: 

m— 1 



E Kj 2 
Kj Xij 

i=0 3 



Yl TF X %^i =0,771-1, 



K 



m—1 



u- x m x m = ^7 x ih x ih^h = 0,m - l,j 2 = 0,n, ji < j 2 , 



i=0 11 



K h 
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Having Kj ^ and K { = 0, all terms, starting with i = m equals to zero. So, matrix A is 
upper orthogonal of size m x m and matrix B is free of size (n — m + 1) x m. For the last 
n — m + 1 columns upper orthogonality has form: 

n V n XV h 

n i 2 \ -< Llp=m ^ 



i=0 - 1 i=m J- lp=m n P 

n U n FT k 

J^i \ ^ LLr>=m 



•v - *.p=m n P 2 w ' 

x a — / ^ x ii ,yj = m,n, 

ll 3 i= m l\.v=m K P 



i=0 Jl j=m 1 lp=m ft P 



because if., = 0. It means the matrix C is upper orthogonal of size (n — m + 1) x (n — m + 1) 
and the matrix O is obligatory zero one of size m x (n — m + 1) (otherwise elements of M 
aren't finite). 

It's easy to verify that inverse matrix has form: 

M ~ {-C-'BA- 1 C- 1 J 

This way of calculating the inverse matrix can easy be generalized to either number of null 
characteristics. 

We will name upper orthogonal matrixes (which also are lower orthogonal) generalized 
orthogonal. We will use the term orthogonal matrix meaning generalized orthogonal matrix 
if isn't stated otherwise. As we can see, the orthogonal matrix set is closed in respect of 
multiplying, it contains the unit element and for any element it contains its inverse. So the 
orthogonal matrix set form isomentry group of space. All motion matrices are generalized 
orthogonal. 



2.5 Orthogonal Matrix as Product of Rotations 

Let X will be orthogonal matrix. We will search the rotation matrices, the product of which 
gives X. Note that The matrix X 9^,(0) have all columns x p of X except % and j ones. These 
columns are x\ = XjCj + i...j(0) + XjSi + \,,,j{<p) and x'j = —^XiSi + i,,.j(4>) + XjCi + i,,,j. 

For the last row let separate elements Xj, % = 0, n in three categories: having characteristics 
K n /Ki equals to 1, and —1. Note that for i row the i element is always of the category 1, 
because its characteristic is Ki/Ki — 1. We will multiply X on the right by *R in ((j)),i = 0, n in 
order to have in the n-th row a single element of category 1 and single element of category — 1, 
different from 0. All these rotations are elliptic ones. For elements of the same characteristic 
x ni and x nj we can use cos0 = , x ni „ and sin0 = , x , nj = . Moreover, always x nn ^ 0. 

J \/x -+x ■ \/x -*-•>• 



X nj 



Now, we have one element of category 1 and one of category —1, different from zero (the 
n-th one and, for example, the p-th one) and element of category 1 has absolute value greater 
then the element of category —1 because for this lower n- normalized row have place equality 
xl„ — xl n — 1. It means that exists d> e R so that cosh <b = , Xnn = and sinh <b = . Xnp = 

and hyperbolic rotation that transforms the element of category — 1, x np in x' np = and the 
element of category 1, x nn in x' nn ^ 0. 
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For category there exist parabolic rotations, that preserves the element of category 1 
(x nn ) and elements of category transform in 0. For this case, if one this element is on q-th 
column, = —x nq . The last non-zero element x nn equals to 1 or —1, because the last row is 
lower n-normalized. 

We can consider the first n columns as having n elements (the last one equals to zero). 
They form orthogonal matrix of size n. The last, (n + l)-th column (without the last element) 
is upper in-orthogonal to first n columns, i = 0, n — 1. As these columns have n elements 
each, (n + l)-th column is obligatory null (excluding the last its element). 

In this stage we can consider the resulting matrix as having size n instead of n + 1 and 
repeat the process for it. Finally obtain the matrix E which has elements on main diagonal 1 or 
— 1 and all rest elements 0. It is the reflection matrix on a point or line or plane or hyperplane. 
Obtain the equality: X T\ q j=1 fflj = E. It's easy to see that X = E nj =g 9JT" 1 (q = (n+l)n/2). 
Strictly speaking, the matrix E can't be presented as product of rotations. In order to identify 
motions of B n with orthogonal matrices, we should name E (which preserve vector product) 
motion. However, these motions are improper (there is no continuous parameterization of 
motion 9Jt(a) on a segment [0,1] such that 9Jl(0) = 7,971(1) = E and all 971(a) are motions 
on all a £ [0, 1]). Having in expression for X determinant of matrices 9Jtj equals to 1 and 
determinant of E is ±1, determinant of X equals ±1. 

2.6 Coordinate and State Matrix 

Consider in some space B ra n+1 vectors v , ...v n . Let coordinates of Vi are [v 0i : ... : v ni ] , i — 0, n. 
Let vectors u$ are ordered and form basis of B n (not obligatory orthonormal). Compose the 
matrix V with elements Vy,i,j = 0, n. Will name it coordinate matrix for vectors v j. Construct 
also the matrix M of size (n + 1) x (n + 1) with elements m^j = Vl ^ V3 . We will name the 
matrix M state matrix of Vi. Having Vi is space basis, elements rriij are all finite. State matrix 
shows how orthonormal is some vector family. It tends to unite one when vectors are more 
normalized and orthogonal to each other 2 . We will demonstrate that volume of parallelepiled 
constructed on vectors Vi equals to | det V\ and 

detM = (detV) 2 . (2.20) 



First, let Vi, i — 0, n are orthonormal. Then the parallelepiped volume is 1, the matrix V is 
orthogonal one and det V = ±1. So, | det V\ equals to parallelepiped volume. All elements on 
main diagonal ma = 1, because all vectors are upper z-normalized. All elements above main 
diagonal are ra^- = 0, i < j, because all vectors u, and Vj are upper r/'-orthogonal (elements 
under the main diagonal may differ from 0). It means that the matrix M is lower triangular 
with all elements on main diagonal equls to 1 and det M = 1 = (det V) 2 . 

Further, note that xQ(y+z) = xQy + xQz, (x+y)Oz = xQz+yQz, (ax)Oy = a(xQy) = 
x (ay),Vx,y, z £ M n ,a £ K. Matrix determinant equals to zero if it contain proportional 
columns or rows. When some row or column of a matrix is multiplied by a, the resulting 
matrix determinant is a times original matrix determinant. When some row or column is sum 



2 If space specification contains null characteristics, some elements of state matrix can have any value, even 
for orthonormal vector family. 
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of two rows / columns, then the matrix determinant equals to sum of determinants of matrices 
containing the first and the second row / column. 

Second, let instead of some Vi use v[ = av^. In this case parallelepiped volume grows in a 
times and | det V'\ = \a\\ det V\. Moreover, det M 1 = a 2 det M — (a det V) 2 = (det V'f . 

Third, let instead of some Vi use v\ — Vi + avj. In this case parallelepiped volume remains 
unchanged, as well as determinant of V, and det M' = det M = (det V) 2 = (det V') 2 . 

Finally, observe, that all matrices V result from orthogonal matrices using operations form 
the second and the third step. It means the equation (2.20) is true for all matrices and volume 
of parallelepiped constructed on vectors Vi equals to | det V\. 

When somebody calculates the parallelepiped volume it's usefull to use the state matrix. 
Its elements don't change on motions and it is always square, even when the number of vectors 
is less then the space dimension (the matrix V isn't square in this case). 

2.7 Plane definition and Specification. Lineals and their 
Specification 

Having B m C IRP m , all m-dimensional planes B n (m < n) lie in IRP n with one global condition 
for vectors X G IRP m C PIR n : X X = 1. Leaving this condition (it doesn't change on any 
motion), we can consider m-dimensional planes of B n as m-dimensional planes of IRP n . 

By definition, m-dimensional plane L m results from subspace B m on some motion. Sub- 
space B m has the first m + 1 columns of unite matrix with dimension n + 1 as its basis. 
Multiplying the basis matrix of B m by some orthogonal matrix result basis matrix of L m as 
first m + 1 columns of orthogonal matrix. Being a subspace, specification of B m contain the 
first m characteristics of specification B n . 

What happens if we take any m + 1 columns of some orthogonal matrix as basis? Let 
column indices i , ix, ■■■i m and i p = 0, n. It's easy to see that motions that preserve this figure 
only change these columns (interior figure motions) or change no these columns (motions of 
B™ that preserve all its points). Thus, figure characteristics K' = K ip ,p = 0,m, or k' p = 

K' v jK' v _ x = Ki v jKi v _ x = rij=ip_i+i kj,p = 1, m is its specification. These figures generally 
speaking are not planes. We will name them lineals. We will name planes also lineals. 

It may happen, that some lineal has K' Q ^ 1 (space and planes have it equals to 1). In this 
case lineal may not intersect the space sphere and may not have image. We will name lineals 
that have no image improper. Although they have no image, their properties help studying 
the space geometry. 

One more interesting case is when the space specification has characteristic —1 and some 
lineal is constructed on limit vectors for this characteristic. Such lineals can't be constructed 
from matrices get as finite product of motions. They can be constructed as limit of infinite 
products. These lineal specifications can't be deduced from space specification. 

For example, let space B 2 has specification {1, —1}. Vectors [0:1:1] and [0 : — 1 : 1] can't 
result from coordinate vectors on finite product of motions. However, there exist translations 
along these vectors (interior lineal translations): 

/ 1 -25 (f) -25 (|) \ 

971= 2S (I) 1-2SH1) -2So 2 (f) , 
V-2S (f) 2SKI) l + 2S 2 (f)/ 
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/ 1 25 (f ) -25 (f) \ 

2B = -25 (*) l-25 2 (f) 2S 2 (f) . 
V-25 (f) -25 2 (f) l + 2S 2 (f)/ 

These motions matrices use functions C (x) and S (x) have characteristic 0, despite the fact 
the space specification doesn't contain zero. These translations are border space motions 
between elliptic and hyperbolic ones. 

2.8 Projection of Vector on Lineal and on its Orthogo- 
nal Completion 

We will name some vector v' projection of vector v on lineal L m , if v' G L m ,v — v' _L L m . Let 
lineal L m is constructed on vectors l , ...l rn . Then v' = Y^T=o v ^h- Evident, v' G L m . Let's 
see: 

m ^ ; 
i=0 * 

= uoZj -— ^ftoy 



for all j = 0,m, in other words, v" — v — v' _L L m . When some i^j = 0, expression has 
undefined value. It happens when some vector direction is orthogonal to all others. In this 
case there is impossible to determine unique orthogonal vector. However, any value of this 
expression, for example 0, is valid, as it corresponds to some orthogonal vector. 



2.9 Basis Change in Lineal. Unique Form of Lineal 

Let L m C B n is some space lineal, defined by matrix of size (n + 1) x (m + 1). The matrix 
columns l^i = 0,m form basis of lineal. Consider vector a — [a : ... : a m ] G L m . Let vector 
coordinates in B m are v — [v : ... : v n \. Then v = L m a. Let 971 be interior motion of lineal 
L m , defined by matrix of size (m + 1) x (m + 1). And let coordinates of a in new basis L' m 
are b = [b : ... : b m ]. Then b = VJla. Now v = L' m b. Having the fact the coordinates of vector 
v in B n don't change, result matrix equality: 

L m a = v = L' m b = L' m {ma) = {L' m m)a. 

This equality doesn't depend on vector a, then 

L rn = L lm m 



is equation of basis change in lineal. 
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It is necessary to find the unique form of lineal definition. Consider the following algorithm 
for the unique basis search: 



1. Let i p ,p = 0, n is basis of B n . Start with empty basis of L m . 

2. Until new basis has less then m + 1 elements, search for i' p as projection of next i p on 
L m . 

(a) If projection isn't null, find new vector i" as projection of i' p on orthogonal comple- 
tion of existing basis U. 

(b) If i" v isn't null, find its position as free index < q < m so that r 2 

i" 

(c) Norm it and add to existing basis l q = —. 



1 i''Oi" > 0. 



2.10 Measure Calculus Between Lineals 




Let X p ,Y q ,p < q < n are two lineals. Let 
Xi,i = 0,p is the basis of X p . Let X' p be 
projection of X p on Y p (Figure 2.3) and let 
x\ be projection of Xi on Y p (they are not or- 
thonormal) . If the volumes of parallelepipeds 
constructed on vectors Xi and x\ are equals 
to V x and V' x respectively and the angle be- 
tween X p and Y q is measurable and equals 
to 0, then has place the equality: 

v' x = v x c{4>). 

This equality is a particular case of (2.8), 
when k\ = 0, T(x) = x. In our case al- 
ways k\ = 0, because the space model is lin- 
ear. As were discussed earlier, V x = 1 and 
V x = ^det M' x , where M' x is state matrix of 
vectors x'y. 



C((f>) = ^detM^. 



(2.22) 



Figure 2.3: Measure calculus between lineals 
X 2 and Y 2 . 



It may happen that characteristic of <fi equals 
to zero and we can't calculate C~ l (<p). If we 
project vectors Xi on Y^ orthogonal comple- 
tion of Y q (suppose it has dimension at least 
p), get X" p constructed on vectors x" with 
the volume V x , then by (2.13) get: 

v; = vjs {</>). 

Or, having M" state matrix for vectors x ; -, 



S((/)) = ^detM'J. 



(2.23) 



2.11. VOLUME CALCULATION 
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If the dimension of Yf[ is less then p, then we can get S((f>) by projecting of Y± on X p . If <f> 
isn't measurable, then the angle if) between X p and Yj_ is measurable and: 

S(if>) = y^tW x , (2.24) 
C(ip) = v/det M». (2.25) 

The angles (f> and ^ present measure between lineals X p and F 9 . Having values C((f>), S((f>) 
and C (if)) , S (if)) , it is possible to determine <f> and '0. The measure characteristic of <f> and 
if> equals to measure characteristic between X' p and X" p . Depending on this characteristic, 
situation can be one of the following: 



If characteristic equals to 1, then detM^. + det M" = 1 and <f> = tan 1 J , if) 



tan" 1 / detM * 



det M" ' 



If characteristic equals to 0, then either det M x = 1 and <f> = ^/det M xl if) = oo, or 
det M" = 1, and = oo, ^ = ^/detM£. 



det M" 



If characteristic equals to —1, then either det M' x — det M" = 1 and <f) = tanh 1 

■0 isn't measurable, or det M x — det M' x = 1 and <f> isn't measurable, = tanh -1 ' ' ' 1 
or det M" = det M' x and <p = ip = oo. 



det M' ' 



det Mil ' 



2.11 Volume Calculation 

We can see that for any B n seen as a unit sphere in IR n+1 , the surface is orthogonal to radius. 
Let I,F 6 I" and the distance between X and Y is small. Let O = (0,0, ...,0) is origin of 
]R n+1 . We will see that (O — X)Q(Y — X) = when Y — > X in sense of distance between them. 
O-X = -X, (O-X)Q(Y-X) = -XQ(Y-X) = XQX-XQY = l-C x (d(X,Y)), where 
d(X, Y) is distance between X and Y. When Y -)• X, d(X, Y) -» and 1 - Ci(d(X, Y)) -> 0. 

Let A, 5 e B 1 . A = [Ci(a) : Si (a)], 5 = [Ci(/3) : Si (/?)], where Ci(a;), S'i(x) and Ti(ar) 
are defined as (1.2), (1.3) and (1.4). Let calculate in M 2 the area of B 1 sector between A and 
B. In Euclidean polar system the argument tan0 e = y/x = Si(<j>) / C x ((f>) = Ti(0), where <f> 
is native argument in B 1 . The Euclidean radius p = \J x 2 + y 2 = a/ Cf((f>) + S 2 ((f>). Having 

^ e = (l+T?{<f>))Cf(<p)' tne area is: 

=2/, omrm ^^L »=2* a =—- 

That is, 25 equals to length AB. 

Let F C B" be some figure with volume (in sense of B n ) Vb- We will name the volume 
(in sense of IR n+1 ) of cone with base F C B n and vertex O ^ B n origin of IR n+1 (Figure 
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ra+l 



2.4). As B n is orthogonal to radius, F also is orthogonal. The radius equals 1, because 
VX e W,X © X = 1. Then for each figure F cM n have place equality: 

V B = (n + 1)V R 

As motions preserve B n and the absolute value of their matrices' determinant is 1, all 
motions preserve Vr and thus, they preserves also Vb. 



Chapter 3 

Theory Application 



3.1 Space and Lineal Specification Search Algorithm 

As we can see, the theory described in this book is universal and easy applicable. However, one 
issue stops somebody from using it. Geometric spaces are classified and defined different from 
the way adopted here. Therefore, in order to not loose the feeling of reality, we will describe an 
algorithm aimed to find specification for some geometric space. The algorithm can be applied 
to any space where have sense notions of points, lines, planes, subspaces, distances, angles and 
/ or motions. 

1. Let m equals to the greatest number of general situated points, or same, the lowest 
number of vertices in a polyhedron of positive volume. 

2. Count space dimension as n — m — 1. 

3. Name points O-dimensional planes and lines 1-dimensional planes. 

4. For % — 1, n do: 

(a) If among (i — l)-dimensional planes there are non-congruent ones, then the space 
definition or space terminology is inconsistent. Theory still can be used, however, in 
order to understand it correctly, it is necessary to modify terminology or to define 
otherwise some space elements (about it later). 

(b) If the measure between {i — l)-dimensional planes is bounded, then hi = 1. 

(c) If the measure between (i — l)-dimensional planes is scalable, then fcj = 0. 

(d) Otherwise, fcj = —1. 

5. Having space dimension n and specification {k±, k n }, use theory. 

The necessity of proper terminology, uniform among all spaces is required by wish to have 
such a theory, that isn't misleading and helps us to study the space structure and to compare 
it with other spaces. Still, under inconsistent theory / terminology we should understand it 
has a contradiction, but failing it to match to theory / terminology that is common today. 
We assume the following here: 
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• All the planes of any dimension are congruent, including points and lines. 



• Theory allows the duality principle of (m— 1) -dimensional planes and (n— m)-dimensional 
ones. 



We should mention that 'common terminology' may change over the time. In order to un- 
derstand what it is consider an example of inconsistent terminology. The Minkowskii space is 
successfully used in physics to describe the theory of relativity. Unfortunately, from geometry 
point of view, it have no proper terminology. The notions of 'space-like lines', 'time-like lines' 
and 'light-like lines' have sense in physics, but not in geometry. Corresponding geometric 
notions are: 'I-st category lines', 'II-nd category lines' and 'III-rd category lines'. No space 
motion maps some line of a category into some line of another category. There is no con- 
tradiction here, but there is an inconsistence. What happens if somebody wants to define a 
space with five categories of lines 1 ? Nobody defines several kinds of points. All points are 
congruent 2 . Why shouldn't be lines all congruent? At the other hand, relative position of 
points may differ. If we name lines only the I-st category of lines, then we should exclude 
II-nd and III-rd category of lines from lines. At the first look, it conflicts with the axiom that 
claims any two points can be connected with a line. But this axiom may have no place in 
other spaces. In contrast, just Euclidean geometry, where all points are connectable, gives 
us an example of parallel lines (that have no common point). Using the duality principle, it 
should exist the notion of non-connectable points (that have no common line). 

It should be mention, that even for somebody feels comfortable using this theory, the 
algorithm described earlier may help to determinate the specification of some exotic lineals (for 
example, of ones defined as limit lineals, which aren't deductable from the space specification). 



3.2 Some Special Spaces 

Many linear spaces are defined using the quadric form of distance d 2 (X, Y) = (Y— X)o(Y— X). 
As for these spaces k± — 0, Kq = 1, Ki — 0, i > 0. In this case, the equality 1 = C(d(X, Y)) = 
X Y = 1 is trivial and can't be used for distance calculation. Consider one more vector 
product — such as (X Yf + k x {X Yf = 1, VX, Y G B n , k = {-1, 0, 1}. This product 



depending on concrete space, there can exist more categories of two-dimensional planes. For further 
dimensions of planes the number of their categories grows. 

2 Thc notion of points on infinity is used in projective geometry. These points are non-congruent with 
others. The terminology is not common in a scope of analytic geometry. 
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is similar to exterior vector product. Change 1 = (X X)(Y Y): 
(X®Y) 2 = — ((XqX)(YqY)-(XqY) 2 ) 

= ((e^*) (e W) - (f>^) (E'vo//, 

^ n rt 

= IT E E ~ x i x iViyj) 

Kl i=0 j=0 

1 n 

1 i<j=0 

1 n 

— ^ KiKjfayj - Xjyi) 2 



i<j=0 



So, 



X0F = 



1 n 

— J] KiKjixiVj-Xjyi) 2 . (3.1) 

1 i<j=0 



Note that from l n = {i 6 MP n | x x = x 2 = 1} result x = 1 or x = -1. As x G B n 
implies — x G B n we can consider xq = 1. We will use operator for distance d(X, Y) 
between points X and Y. Note that having C 2 (x) + /c 1 S' 2 (X) = 1, Vx G R, fci = { — 1, 0, 1} and 
(X0X) 2 + A; 1 (X0F) 2 = 1, VX, F G B n , A* = {-1,0,1} result S(d(X,Y)) = X ®Y,VX,Y G 
B n : 

1 n 

d 2 (X,Y) = S 2 (d(X,Y)) = (X ®Y) 2 = - KiKjixiVj - x^ 2 

1 i<i=0 

In this sum all non-zero terms are those for which % — 0: 

d 2 (X,F) = l^^(x 0% -^yo) 2 

i=i 1 

n n 
j=l p=2 

In this equality don't appear x or y . We can consider B n a hyperplane of IR™ +1 with equation 
x = 1 and specification {k 2 , ...k n }. We can identify it with W 1 . Then the equality above is 
equivalent to (Y - X) © (Y - X). 

It means firstly, that scalar product of vestors in lnear spaces (ki = 0) induces the same 
metrics that is used in the model, and secondly, that non-linear spaces with specification 
{ki, k 2 , ...kn} (ki 7^ 0) are best approximated by linear spaces with specification {0, k 2 , ■■■k n }. 
Note also that from here deduce that non-linear space with specification {ki : ...k n } is enclosed 
in model meta-space of greater by one dimension, of which specification is {0, hi, ...k n }. 
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We can use this quadric form to search for all characteristics except k\ = 0. We will use 
this method in order to describe some special spaces by specifying their specifications. 

Case 1. Elliptic, Euclidean and Hyperbolic Spaces. Elliptic, linear (Euclidean) and 
hyperbolic (Bolyai-Lobachevsky) spaces have characteristic k\ equals to sign of space curvature 
ki — 1 for elliptic space, k± — for linear space and k = — 1 for hyperbolic space. 

All these spaces are usually approximated by Euclidean one. We can calculate the rest of 
characteristics using euclidean quadric form. Let dimension is 3: 

= (yi - xxf + k 2 (y 2 - x 2 ) 2 + k 2 k 3 (y 3 - x 3 ) 2 
so k 2 = 1 and k 2 k 3 — 1, k 3 — 1. 

Case 2. Minkowskii Space. The distance between X and Y is calculated (for time-like 
vectors) as 

d 2 (X,Y) = (yi -xif - (y 2 - x 2 ) 2 - (y 3 - x 3 ) 2 - (y A - x 4 ) 2 

= {yi - xi) 2 + k 2 (y 2 - x 2 ) 2 + k 2 k 3 (y 3 - x 3 ) 2 + k 2 k 3 k 4 (y 4 - x 4 ) 2 

where coordinate 1 is time-like and coordinates 2, 3 and 4 are space-like. So k 2 = — 1, 
k 2 k 3 = —1, k 3 — 1 and k 2 k 3 k 4 = —1, k 4 = 1. As Minkowskii space is linear, k\ = 0. 

If we introduce curvature in space, its structure changes. For example, let ki — 1. Then 

X Y = x y + hxiyi + k x k 2 x 2 y 2 + hk 2 k 3 x 3 y 3 + k^k^x^ 
= x y + xiyi - x 2 y 2 - x 3 y 3 - x A y A 

so time characteristic becomes elliptic and space characteristic becomes hyperbolic. If k± — 1, 
then 

X Y = x y + fciXiyi + k x k 2 x 2 y 2 + k x k 2 k 3 x 3 y 3 + k x k 2 k 3 k A x A y A 
= x y - xiyi + x 2 y 2 + x 3 y 3 + x A y A 

and time characteristic becomes hyperbolic and space characteristic becomes elliptic. 

Case 3. Minkowskii Space with 2-dimensional Time. Consider a 4-dimensional space 
with distance quadric form that has 2 positive signs and 2 negative. This space is sometimes 
named Minkowskii space with 2-dimensional time: 

d 2 (X,Y) = { yi - Xl ) 2 + {y 2 - X2 f-{y z - Xz f-{ yi - Xi f 

= (yi - xi) 2 + k 2 (y 2 - x 2 ) 2 + k 2 k 3 (y 3 - x 3 ) 2 + k 2 k 3 k A (y A - x A ) 2 

So k 2 = 1, k 2 k 3 = —1, k 3 = —1 and k 2 k 3 k A = —1, k 4 = 1. As for all lineal spaces, ki — for 
it. 
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Case 4. Spaces with Degenerate Distance Quadric Form. Consider linear 4-dimensional 
space (ki = 0) with degenerate distance quadric form: 

d 2 (X,Y) = ( Vl - Xl f + (y 2 - x 2 ) 2 + (y 3 - x 3 ) 2 

= (y 1 - xi) 2 + k 2 (y 2 - x 2 f + k 2 k 3 (y 3 - x 3 ) 2 + k 2 k 3 k 4 (y 4 - x 4 ) 2 

so k 2 = k 3 = 1 and k 4 = 0. 
It means that motions: 

x[ = Xi 

x 2 = x 2 

x' 3 = x 3 

x 4 = (f)iX 1 + (f) 2 X 2 + <f) 3 X 3 + X A 

are all valid. 

However, transformation: 

x\ = Xi 

x 2 = x 2 
x' 3 = x 3 

x' A = (j) 4 X 4 

is not a motion. Although it preserved distance, it doesn't preserve volume except <f> 4 = 1 or 
— 1. It is an example of angle scaling. 

3.3 Spaces as Product of their Subspaces 

Another way to define spaces is by product of their subspaces. It is necessary to be accurate 
here. The geometric space isn't only a structure of points. It is also the structure of all its 
subspaces. It is mistake to think that having IR 1 is isomorphic to one-dimensional Euclidean 
space E 1 , from R 1 x R 1 = R 2 results E'xE 1 = E 2 (using specification notation, {0} x {0} = {0, 
1}). The problem is the product doesn't define way to measure the angle between multiplied 
subspaces. It can be defined in several ways, for example, {0,0} or {0, —1}. 

The situation is even worse when multiplied subspaces with different specification X m and 
Y n . One-dimensional images can be constructed in two ways: X 1 x Y° (isomorphic to X 1 ) and 
X° x Y 1 (isomorphic to Y l ). And if X 1 and Y l have different specifications, these two one- 
dimensional lines aren't congruent. For example, if somebody wants to construct geometry on 
a cylinder, first thing he or she thinks of is S 1 x E 1 ({1} x {0}), where S 1 is one-dimensional 
elliptic space. In this case some lines are circles, some are lines and others are right and left 
helices, that may not intersect, intersect in one point or intersect in infinity of points. As an 
example of complete geometry on a cylinder you can take the space with specification {1,0}. 

Additionally, one should not consider that if from algebraic point of view E 1 is isomorphic 
to H 1 (one-dimensional hyperbolic space), then constructions like H 2 x E 1 ({ — 1, 1} x {0}) 
and H 2 x H 1 ({ — 1,1} x { — 1}) are also isomorphic. From geometric point of view, E 1 is 
scalable, while H 1 is not (the mutual departure of points is possible, however it can't be linear). 
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In contrast, it is possible to construct spaces with specifications { — 1,1,0} and { — 1,1,-1}, 
which differ one from another by the fact that in first one on a twodimensional plane doesn't 
containing some point there is the only point that isn't connectable with it, and for the second 
space the number of such a points is infinity. 
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